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I.  INTRODUCTION 


The  electromagnetic  properties  of  the  Earth’s  ionosphere  are  important  to  many 
military  applications,  such  as  HF  radio  communications  and  over-the-horizon  radar.  These 
properties  depend  on  the  ionosphere’s  electron  density  profile,  which  varies  over  time, 
altitude,  and  geographic  location.  Currently,  ground  based  ionosondes  arc  used  to  measure 
the  ionospheric  electron  density.  This  technique  is  accurate  and  fast,  but  spatially  limited. 
Global  electron  density  maps  arc  based  on  the  profiles  from  only  a  few  sites.  True  global 
measurements  are  desired.  As  satellite-based  ionosondes  are  impractical,  efforts  arc  being 
directed  toward:  other  potential  worldwide  sensing  methods. 

Observation  of  the  Earth’s  airglow  spectra  is  such  a  technique.  The  principal  ions 
in  the  lower  thermosphere  are  chemically  coupled  to  the  various  neutral  molecules  and 
atoms  present.  The  densities  of  the  neutral  species  can  be  determined  from  measurements 
of  the  ultraviolet  emission  spectrum  of  the  ionosphere.  Electron  densities  can  then  be 
inferred  from  these  measurements. 

This  thesis  presents  an  analysis  of  the  ultraviolet  emissions  of  nitric  oxide  (NO)  in 
the  Earth’s  atmosphere.  These  emissions  were  observed  by  the  Naval  Postgraduate 
School’s  Middle  Ultraviolet  Spectrograph  (MUSTANG)  experiment  flown  aboard  a 
NASA  sounding  rocket.  The  experiment  is  part  of  an  ongoing  effort  to  develop  passive 
remote  sensing  techniques. 
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A.  THESIS  OBJECTIVES 


This  work  focuses  on  thirteen  averaged  spectra  obtained  by  the  MUSTANG 
experiment  over  an  altimde  range  of  100  to  200  kilometers.  Nitric  oxide  emissions  with 
wavelengths  from  2000  A  to  2500  A  are  examined.  Each  observed  spectrum  is  modeled 
by  a  synthetic  spectrum  in  which  the  relative  strengths  of  the  most  prominent  gamma  and 
epsilon  band  transitions  are  varied  independently.  Eleven  emission  features  of  NO  are 
identified  and  analyzed  in  the  wavelength  range  of  interest. 

The  goal  of  the  analysis  is  to  determine  the  column  density  of  nitric  oxide  at  each 
altitude.  A  first  step  in  this  process  is  to  obtain  a  temperature  profile  from  the  data.  These 
temperatures  are  used  for  the  rest  of  the  fitting  procediue. 

B.  THESIS  OUTLINE 

The  thesis  is  divided  into  six  chapters.  Chapter  n  gives  some  background  on  the 
stratification  of  the  atmosphere  as  well  as  some  fundamentals  of  ni*Tic  oxide  chemistiy- 
in  the  lower  ionosphere.  In  Chapter  HI,  a  brief  review  of  diatomic  molecular  spectra  is 
presented.  Application  of  the  results  to  the  calculation  of  the  synthetic  spectra  is 
discussed. 

Chaf  .er  IV  describes  the  experiment  and  the  MUSTANG  instrument  used  in 
obtaining  the  spectra.  The  data  are  also  presented.  These  data  are  analyzed  in  Chapter  V. 
Each  observed  spectrum  is  fitted  with  a  synthetic  spectrum  as  closely  as  possible.  Causes 
for  variations  between  the  model  and  the  observed  spectra  are  discussed.  Chapter  VI 
concludes  with  a  summary  of  the  findings  and  suggestions  for  further  research. 


2 


n.  BACKGROUND 


A.  INTRODUCTIOi> 

Although  nitric  oxide  is  a  trace  element  in  the  atmosphere,  it  plays  a  major  role  in 
atmospheric  chemistry.  In  the  stratosphere  and  mesosphere,  it  has  ’^en  noted  that  NO  acts 
as  a  catalyst  in  the  dertruction  of  ozone  (Oj).  The  low  first  ionization  poteniial  of  nitric 
oxide  makes  it  an  important  source  of  ions  in  u.e  lower  ionosphere.  The  ionization  of 
nitric  oxide  by  the  solar  Lyman-a  mission  at  1216  A  is  responsible  for  the  formation  of 
the  D  region  of  the  ionosphere.  A  brief  description  of  the  earth’s  atmosphere,  ionosphere 
and  nitric  oxide  chemistry’  follows. 

B.  THE  NEUTRAL  ATMOSPHERE 

Variations  in  the  Earth’s  atmosphere  are  caused  by  differential  heating  from  solar 
radiation,  convective  mixing,  and  the  rotation  of  the  Earth.  There  results  a  limited 
stratification  of  the  atmosphere  which  is  commonly  described  in  two  ways;  either  in  terms 
of  the  temperature  gradients  in  the  atmosphere  or  by  the  chemical  composition. 

As  shown  in  Figure  2-1,  temperature  gradients  divide  the  atmosphere  into  four  basic 
regions.  The  troposhere  occupies  the  lowest  ten  kilometers  of  the  atmosphere  and  is 
characterized  by  a  steady  decrease  in  temperature  with  height.  It  is  in  this  region  that  our 
"weather"  occurs.  After  passing  through  an  isothermal  region,  the  temperature  begins  to 
increase  with  altitude  up  to  45  kilometers.  This  region  is  known  as  the  stratosphere. 
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Above  this,  in  the  mesosphere,  the  temperature  again  decreases  to  about  180  K  at  an 
altitude  of  5  kilometers.  This  is  the  coldest  region  of  the  atmosphere.  Above  the 
mesosphere,  the  kinetic  temperature  of  the  atmosphere  increases  with  altitude  towards  a 
value  dependent  on  solar  activity.  This  region  is  kno.vn  as  the  thennosphere 


KINETIC  TEMPERATURE  l°KI 

Figure  2-1  Temperature  and  Composition  of  the  Atmosphere  (Tascione,  1988) 

The  lower  boundary  of  the  thermosphere  marks  a  division  between  two  chemical 
composition  regions  of  the  atmosphere.  The  troposhere,  stratosphere,  and  mesosphere  all 
make  up  the  region  known  as  tlie  homosphere.  The  homosphere  is  composed  of 
approximately  78%  nitrogen,  21%  oxygen,  and  1%  of  other  gases.  Mixing  in  the 
homosphere  keeps  the  relative  densities  of  these  gases  nearly  constant  throughout  the 
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region.  Above  100  kilometers,  the  combination  of  low  densities  and  the  stable  temperature 
gradient  cause  mixing  to  cease.  Each  atmospheric  contituent  reaches  its  own  separate 
hydrostatic  equilibrium.  The  relative  densities  of  the  gases  change  with  altitude,  giving 
this  region  the  name  heterosphere.  From  500  kilometers  to  an  indefinite  Limit  lies  the 
exosphere,  a  region  of  high  kinetic  energies  and  low  densities. 

C.  THE  IONOSPHERE 

Above  50  kilometers,  photoionization  and  low  recombination  rates  produce  a 
partially  ionized  region  known  as  the  ionosphere.  This  region  is  divided  on  the  basis  of 
free  electron  densities  into  the  D,  E,  F,,  and  Fj  layers.  These  layers  form  due  to  changes 
in  the  composition  and  solar  absorption  characteristics  of  the  atmosphere.  Electron 
concentrations  in  the  ionosphere  vary  with  solar  activity.  The  ultraviolet  and  x-ray  output 
of  the  Sun  has  been  conelated  with  the  solar  flux  at  2800  MHz.  Commonly  called  FI  0.7, 
this  value  is  measured  daily  and  is  taken  as  a  standard  indicator  of  solar  activity. 

The  D  region  of  the  ionosphere  extends  from  about  50  to  90  km  in  altitude.  As  the 
lowest  lying  region  of  the  ionosphere,  it  is  produced  by  the  deepest  penetrating  ionizing 
radiation  incident  on  the  atmosphere.  The  major  source  of  ions  is  from  the  ionization  of 
NO  by  the  Lyman-a  emission  at  1216  A. 

The  E  region  is  generally  taken  to  lie  between  90  and  130  km.  The  electron  density 
is  due  primarily  to  the  photoionization  of  molecular  oxygen  by  Lyman-P  and  and  soft 
solar  x-rays.  It  has  a  daytime  electron  density  two  orders  of  magnitude  higher  than  the 
D  region’s. 
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In  the  F  regions,  the  principal  ion  is  O*,  formed  by  the  photoionization  of  atomic 
oxygen.  Molecular  oxygen  and  nitrf'gen  also  contribute.  The  region  is  characterized  by 
a  relatively  constant  electron  density.  In  the  Fj  region,  which  lies  above  160  km,  the 
electron  density  reaches  its  maximum  value. 

The  presence  of  nitric  oxide  in  the  lower  ionosphere  is  a  result  of  the  interaction 
of  solar  radiation  with  the  atmosphere.  Processes  which  lead  to  the  dissociation  of  Nj  are 
the  source  of  what  is  known  as  the  odd  nitrogen  family.  This  family  consists  of  all 
photochemicaUy  active  nitrogen  containing  species.  Of  these,  nitric  oxide  is  the  most 
abundant. 

The  mechanisms  leading  to  the  production  and  loss  of  odd  nitrogen  will  not  be 
explored  in  this  thesis  ,  but  the  chemical  reactions  have  been  tabulated  by  Cleary  (1986). 
Figure  2-2  is  a  schematic  representation  of  these  reactions  in  the  thermosphere.  The  figure 
sho  ws  that  the  production  of  odd  nitrogen  is  a  photolytically  driven  process.  Photons  and 
photoelectrons  interact  with  the  major  constituents  of  the  atmosphere,  forming  O/,  O"^, 
and  N/.  In  turn,  chemical  interaction  of  these  ions  with  the  atmosphere  and  the  solar  flux 
produce  the  odd  nitrogen  family.  Many  of  these  reactions  terminate  at  NO"^,  whose  only 
significant  loss  mechanism  is  dissociatve  recombination.  Figure  2-2  also  shows  that  while 
N(^D)  is  a  source  for  neutral  NO,  N(^S)  is  both  a  source  and  a  sink.  The  abundance  of 
nitric  oxide  is  sensitive  to  the  production  rate  of  N(^D). 
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IONOSPHERIC  ODD-NITROGEN  PHOTOCHEMISTRY 


Figure  2-2  Ionospheric  Odd  Nitrogen  Photochemistry 


m.  DIATOMIC  MOLECULAR  SPECTRA 


A.  INTRODUCTION 

Spectroscopic  analysis  is  a  useful  tool  which  has  been  applied  to  a  variety  of  fields. 
From  the  spectral  emission  and  absorption  features  of  a  gas,  one  can  determine  its 
chemical  composition  and  physical  conditions.  Study  of  the  Earth’s  atmospheric  dayglow 
provides  a  means  to  map  the  ionosphere  in  terms  of  composition,  temperature,  and 
reaction  processes  for  various  altitudes  and  geographic  latitudes. 

To  better  understand  the  data  and  the  development  of  the  synthetic  spectra,  a  brief 
review  of  diatomic  molecular  spectra  is  presented.  A  very  thorough  treatment  of  this 
subject  is  given  by  Herzberg  (1950). 

B.  DIATOMIC  MOLECULES 

The  simplest  model  of  the  diatomic  molecule  is  a  dumbbell  in  which  two  spheres 
of  mass  M,  and  Mj  are  separated  by  a  distance  and  bound  by  electrostatic  attraction.  In 
addition  to  the  electronic  motion,  the  molecule  has  two  possible  modes;  vibrational  and 
rotational  motion.  From  quantum  mechanics,  it  is  knowm  that  these  motions  cannot  take 
on  arbitrary  values,  but  rather  are  quantized.  The  energies  of  the  motion  can  be  expressed 
as  term  values,  where  the  temi  value  is  the  energy  divided  by  he.  The  rotational  term 
value,  F(J),  is  given  approximately  by  F(J)  =  BJ(J+1),  where  J  is  the  rotational  quantum 
number  (J=0,l,2,...)  and  B  is  a  constant  of  the  molecule.  Likewise,  the  vibrational  term 


8 


value  is  given  approximately  by  G(v)  =  Q)(v4^),  where  v  is  the  vibrational  quantum 
number.  These  approximations  do  not  account  for  the  anharmonic  or  non-rigid  nature  of 
the  actual  motion.  However,  for  a  given  electronic  state  with  vibrational  level  v  and 
rotational  level  J,  the  term  value  for  the  combined  motion  can  still  be  written  as 

r  =  7,  G(v)  +F(J)  (3-1), 

where  is  the  electronic  energy.  The  electronic,  vibrational,  and  rotational  states  of  a 
molecule  combine  in  a  nested  manner.  Figure  3-1  represents  the  relationship  between  the 
energy  levels. 

The  orbital  configmation  of  the  molecule’s  electrons  determine  its  electronic  energy. 
As  in  atoms,  the  electrons  form  a  total  resultant  orbital  angular  momentum  L  and  a 
resultant  spin  angular  momentum  S.  Since  the  electrostatic  field  of  the  molecule  is 
symmetric  along  the  intemuclear  axis,  L  will  precess  about  the  axis  of  symmetry. 
However,  the  value  of  L  along  the  intemuclear  axis  is  a  constant  of  the  motion  and  is 
denoted  by  (=  L,  L-l,...,-L).  The  magnitude  of  Ml  is  given  the  symbol  A. 

The  resultant  spin  S  is  not  affected  by  the  molecular  electrostatic  field.  It  precesses 
about  the  intemuclear  axis  only  if  there  is  a  magnetic  field  present.  For  A  ^  0,  there  is 
an  internal  magnetic  field,  and  S  precesses  about  the  axis  of  symmetry  with  a  constant 
component,  M,.  Designated  Z,  it  can  take  the  values  S,  S-1,...,  -S,  giving  a  multiplicity 
of  (2S+1). 

The  spectroscopic  notation  for  molecules  is  completely  analogous  to  that  of  atoms. 
The  total  orbital  angular  momentum  is  given  a  Greek  letter  Z,  11,  A,...  corresponding  to 
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Figure  3-1  Electronic,  Vibrational,  and  Rotational  Energy  Levels  (Eisberg  and  Resnick 
1985) 
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L  =  0,  1,  2,..  respectively.  The  multiplicity  of  the  state  is  written  as  a  left  superscript. 
When  Z  is  defined  (A  0),  tlie  value  (A  +  Z)  is  written  as  a  right  subscript.  Furthermore, 
the  electronic  states  are  labelled  by  letters.  The  ground  state  is  traditionally  given  the 
letter  X,  and  excited  states  of  the  same  multiplicity  are  designated  by  A,  B,  C,  etc.  States 
with  diffeient  multiplicity  are  labeled  by  lower  case  letters. 

Figure  3-2  is  a  potential  energy  diagram  for  nitric  oxide  labelled  with  the 
spectroscopic  notation.  The  transitions  of  interest  to  this  analysis  are  the  gamma  and 
epsilon  bands.  The  gamma  band  refers  to  the  transition  between  the  A^Z  and  X^Il  states. 
Similarly,  the  epsilon  band  refers  to  the  EFZ  to  X^H  transition. 

C.  LINE  POSITIONS 

The  line  position,  a,  or  wave  number  of  a  spectral  emission  line  is  given  by  the 
difference  in  term  values  between  the  initial  and  final  states.  The  line  position  of  a 
transition  from  the  (T/,v',J')  state  to  the  (T/',v"J")  state  is  given  by 

o  =  +  G'(vVg"(v'0  + 

where  the  single  prime  refers  to  the  upper  state  and  the  double  prime  refers  to  the  lower 
state. 

The  set  of  transitions  for  which  T/  and  T/'  are  fixed  are  called  a  molecular  band 
system.  Within  a  given  band  system,  the  transitions  for  which  v'  is  fixed  and  v"  varies 
are  known  as  a  v"  band  progression. 
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INTERNUCLEAR  DISTANCE  (A) 


Figure  3-2  Potential  Energy  Diagram  of  NO  (from  Gilmore,  1965) 
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To  arrive  at  the  correct  expressions  for  G(v)  and  F(J),  one  must  account  for  the 
anharmonicity  of  the  oscillations  and  the  mutual  influence  of  the  rotational  and  electronic 
angular  momenta.  The  electronic  and  rotational  angular  momenta  of  the  molecule  can 
couple  in  many  ways.  The  form  of  F(J)  depends  on  the  type  of  coupling.  That 
development  will  not  be  presented  here. 


D.  EMISSION  RATES 

The  volume  emission  rate  for  a  single  molecular  line  can  be  written  as 


E(n'v'j',  n"v"j')  =  N{n'v'j)  A{n'v'j'.  n"v"j")  (3-3). 

where  E  is  the  volume  emission  rate,  N  is  the  fractional  population  of  the  state  defined 
by  the  quantum  numbers  n',  v',  J',  and  A  is  the  Einstein  coefficient  for  spontaneous 
emission  between  the  given  states. 

The  number  of  molecules  in  the  level  defined  by  n'  and  v'  is  called  N(n'v').  In  the 
steady  state,  this  number  is  given  by 


Nin'y^ 


P{n'v) 


(3-4), 


where  P(n'v')  is  the  excitation  rate  to  the  (n'v')  state,  and  the  denominator  is  the  sum 
over  aU  the  loss  processes.  Using  this,  equation  (3-3)  can  be  rewritten  as 


EinW/i'W)  =  Pin'v^ 


(3-5). 


NjnVj)  Afj" 

Nin'v)  53i4^v' 

y" 

The  principal  excitation  mechanism  for  NO  is  fluorescent  scattering.  The  production 
rate  is  proportional  to  the  solar  flux  and  the  effective  cross  section  for  transitions.  The 
branching  ratio  into  the  vibrational  levels  is  given  by  the  Franck-Condon  factors. 
Assuming  that  the  absorbing  molecules  are  all  in  the  ground  state,  this  gives 

P(/iV0  =  nF  —  'yov' 

me 

where  tcF  is  the  solar  flux,  ^  is  the  wavelength  of  the  incident  radiation,  f„  is  the 
oscillator  strength  for  the  transitions  from  the  ground  state  to  n',  and  is  the  Franck- 
Condon  factor  for  the  v=0  ground  state  to  v'  level. 

The  fractional  population  of  a  molecular  state,  N(n'vT)/N(n'v').  is  determined  by 
the  excitation  process.  Calculations  of  this  ratio,  applicable  to  NO  spectra,  may  be  found 
in  Tatmn  (1967). 

Finally,  the  transition  probability  Aj^j»  can  be  related  to  using  the  Honl-London 
factors  SfJ'J").  Then,  A,.j-  is  expressed  as 


A 


•j'j" 


-  ^y'y" 


2(J^1) 


(3-7), 


wher^  2(J'-(-l)  is  a  normalization  factor.  The  Honl  London  factors  give  the  probability  of 
a  transition  from  an  upper  state  J'  to  a  lower  state  J". 
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With 


“vV"  = 


Equation  (3-5)  can  be  rewritten  as 


(3-8), 


N(n'v')  ""  2J^l 


(3-9). 


The  single  scattering  albedo,  Wy-v'.  c^n  be  determined  from  the  Franck-Condon  factors 


(Barth,  1965). 

E.  SYNTHETIC  SPECTRA 


The  elements  described  in  this  chapter  can  be  used  to  generate  synthetic  spectra. 
The  calculation  for  a  molecular  band  requires  that  the  emission  rates  must  be  determined 
for  each  rotational  line  in  the  band.  The  spectral  region  is  divide  into  wavelength  bins. 
As  the  lines  are  calculated,  they  are  added  to  the  appropriate  bin.  This  process  is  repeated 
until  all  the  contributions  have  been  added.  Algorithms  for  this  were  developed  by  Qeary 
(1985)  and  were  used  in  the  analysis  described  below. 
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rv.  THE  EXPERIMENT 


A.  INTRODUCTION 

The  experiment  was  launched  from  White  Sands  Missile  Range  at  1700  GMT  on 
March  30,  1990  aboard  a  NASA  Terrier  Black  Brant  sounding  rocket.  The  payload 
consisted  of  two  instniments.  One,  a  Naval  Research  Laboratory  (NRL)  spectrograph,  was 
designed  to  observe  emissions  in  the  Fj  region  of  the  ionosphere.  The  second 
instrument,  the  Naval  Postgraduate  School’s  (NPS)  Middle  Ultraviolet  Spectrograph 
(MUSTANG),  was  designed  to  measure  the  emissions  jfrom  NO,  Nj,  and  N*  in  the  D  and 
E  regions  of  the  ionosphere.  The  MUSTANG  experiment  produced  the  data  which  is  the 
subject  of  this  thesis. 

The  sounding  rocket  reached  an  apogee  of  320  kilometers,  with  altitude  infonnation 
provided  by  ground  based  radar  tracking  stations.  During  the  flight,  the  rocket’s  attitude 
control  system  (ACS)  maneuvered  the  payload  to  optimize  the  conditions  for  data 
collection  for  each  experiment.  For  the  MUSTANG  experiment,  this  involved  positioning 
the  payload  to  view  the  horizon  with  a  90  degree  observation-zenith  angle  (horizontal 
viewing)  and  a  90  degree  solar-observation  angle.  Approaching  apogee,  the  payload  was 
maneuvered  through  270  degrees  so  that  it  was  pointing  at  the  zerrith  for  part  of  the 
downward  portion  of  the  flight.  At  approximately  200  kilometers,  the  payload  was  again 
rotated  to  point  towards  the  horizon.  This  positioning  resulted  in  two  sets  of  useable  data 
from  MUSTANG.  On  the  up  leg  of  the  flight,  spectra  were  recorded  in  the  altitude  range 
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of  150  .o  200  kilometers.  On  Ote  do»n  leg.  spectra  were  measured  &om  180  to  100 
kilometers. 

B.  THE  INSTRUMENT 

The  MUSTANG  instrument  is  a  1/*  meter  Ebert-Fastie  spectrograph  with  an  off 
aais  telescope.  The  optical  design  is  shown  scheroaticaUy  in  Figure  4-1.  The  entrance  slit 


Figure  4-1  Optical  Design  of  Spectrograph  (from  Samson,  1967) 
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was  5  millimeters  in  height  by  140  microns  wide  producing  a  field  of  view  of  2.3  degrees 
in  tlie  vertical  and  0.006  degrees  horizontal.  The  grating  ruTng  was  1200  lines  per 
millimeter.  It’s  positiori  was  adjusted  to  produce  a  wavelength  coverage  of  1800  A  to 
3400  A.  The  detector  system  consists  of  an  ITT  F4145  microchannel  plate  image 
intensitier  fiber  optically  coupled  to  a  Hamamatsu  S2300-512Q  512  element  linear 
photodiode  array. 

The  image  intensifier  is  composed  of  a  photocathode,  two  microcharmel  plates  in 
cascade,  and  an  output  phosphor  screen.  The  phoiocathode  was  placed  in  the  focal  plane 
of  the  spectrograph.  Light  fi'om  the  mirror  strikes  the  photocathode,  causing  an  electron 
to  be  ejected.  Electrical  potentials  across  the  image  intensifier  accelerate  the  electron 
through  the  microchannel  plates  (MCP).  Figure  4-2  shows  a  block  diagram  of  the  device. 
The  microchannel  plates  are  similar  to  photomultiplier  tubes  in  that  15,000  electrons  are 
produced  as  output  for  eveiy  one  input.  These  electrons  exit  the  MCP  and  strike  the 
phosphor,  producing  photons.  The  image  intensifier  fiber  optic  output  window  is  in 
contact  with  the  fiber  optic  input  window  of  the  photodiode  anay,  resulting  in  negligible 
loss  of  photons  The  output  of  the  photodiode  is  an  analog  voltage  proportional  to  the 
incident  intensity. 

The  wavelength  calibration  or  the  spectrograph  was  determined  by  identifying 
atomic  lines  from  the  spectra  of  platinum  and  mercury.  Th^  sensitivity  of  the  instrument 
was  determined  by  illuminating  a  barium  sulfate  (BaSO^)  screen  using  a  standard 
deuterium  lamp  and  a  tungsten  filament  quartz-halogen  lair.p.  At  the  time  of  this  paper, 
the  post  flight  calibration  was  not  complete.  As  a  result,  the  sensitivity  calibration  is  still 
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Figure  4-2  Block  Diagram  of  the  Image  Intensifier 

preliminary.  However,  the  systematic  error  introduced  by  the  uncertainty  in  determining 
the  sensitivity  is  typicaUy  20%  for  these  types  of  experiments  (Qeary  and  Barth, 1987; 
McCoy, 1983).  The  detector  response  to  an  input  was  found  to  be  Linear  up  to  85%  of  its 
maximum  output.  Wavelength  resolution  of  the  instrument  was  found  by  determining  the 
instrument  slit  function.  This  is  simply  the  instrument  response  to  a  delta  function  input. 
The  delta  function  was  approximated  by  the  atomic  lines  of  mercury.  Wavelength 
resolution  was  found  to  be  10.6  A. 

C.  THE  OBSERVED  SPECTRA 


The  observations  are  of  the  atmospheric  dayglow  under  mid-day,  mid-latitude 
conditions.  The  solar  zenith  angle  was  62.9  degrees  and  the  F10.7  was  186.9. 
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Data  integration  time  was  50  milliseconds  (ms).  Output  from  the  photodiode 
detector  was  converted  into  a  stream  of  512  ten  bit  digital  words.  These  data  were  sent 
via  telemtry  to  ground  stations  where  they  were  recorded.  The  telemetry  signal  also 
carried  mission  elapsed  time  data,  making  it  possible  to  correlate  the  spectral  data  with 
altitude  information. 

The  spectra  were  averaged  in  ten  kilometer  bins.  This  resulted  in  an  average  of  100 
spectra  per  altitude  bin.  The  actual  number  of  spectra  in  each  bin  varied  by  as  much  as 
20%  at  the  altitude  extremes,  due  primarily  to  changes  in  the  vertical  velocity. 

A  total  of  13  average  spectra  were  obtained.  Figures  4-3  through  4-5  are 
reprensentative  of  the  high,  mid,  and  low  altitude  observations.  Qearly  seen  are  the  nitric 
oxide  emissions,  particularly  the  7(1,0)  transition  at  2150  A.  The  analysis  of  the  following 
chapter  concentrates  on  the  NO  emissions  in  the  2000  A  to  2500  A  range. 


20 


300 


O  O'" 


c  c 

N 


(v/a)  ‘4ISU91UI 


Figure  4-3  Spectrum  at  195  km  (Up  Leg) 
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Figure  4-4  Spectrum  at  155  km  (Down  Leg) 
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Figure  4-5  Spectrum  at  105  km  (Down  Leg) 
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V.  DATA  ANALYSIS 


A.  INTRODUCTION 

In  this  chapter,  the  data  are  analyzed  by  comparison  with  synthetically  generated 
spectra.  Eleven  major  features  of  the  observed  data  are  identified  and  fitted  to  the 
synthetically  generated  spectra  by  independently  varying  the  corresponding  synthetic  nitric 
oxide  emissions  in  the  2000  A  to  2500  A  range.  The  objective  of  the  analysis  is  to  extract 
NO  line-of-sight  column  densities  from  the  spectra.  In  addition,  a  temperature  profile  is 
obtained.  Variations  of  the  data  from  model  predictions  are  explored  and  discussed. 

The  data  readout  fi-om  the  instrument  was  dependent  on  telemetry  clocks.  During 
preflight  preparations,  it  was  discovered  that  the  pulse  width  of  one  of  the  NASA  timing 
signals  did  not  meet  specifications.  This  caused  every  seventeenth  data  word  to  be 
dropped  from  the  telemetry  stream.  This  was  corrected  during  the  data  reduction  by 
inserting  an  averaged  value  where  the  missing  data  should  have  been  and  shifting  the 
adjacent  subset  of  observed  data  points.  An  additional  wavelength  discrepancy  was  found 
when  fitting  the  data  with  the  synthetic  spectra.  The  data  appeared  to  be  shifted 
progressively  with  increasing  wavelengths.  Figure  5-1  shows  the  wavelength  shift  as  a 
function  of  band  head  wavelength.  The  low  altitude  spectra  were  used  to  determine  the 
size  of  the  shift.  The  postflight  calibration  is  still  in  progress.  As  a  result,  it  cannot  be 
determined  at  this  time  whether  this  shift  is  an  instrumental  artifact  or  an  actual 
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Wavelength  Shift  (X) 


disagreement  between  theory  and  data.  During  the  fitting  process,  the  data  were  shifted 
by  the  necessary  amount. 


■Wavelength  (X) 

Figure  5-1  Wavelength  Shift  Deduced  From  the  Fit  to  the  NO  Bands 

B.  APPLICATION  OF  THE  SYNTHETIC  SPECTRA 

The  physical  dimension  of  the  synthetic  spectra  is  an  emission  rate  per  atom  per 
wavelength  bin.  The  intensity  of  the  observed  data  is  in  Rayleighs  per  Angstrom  (one 
Rayleigh  equals  10®  photons  cm‘^  sec'*).  Thus,  by  fitting  the  data  with  the  synthetic 
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spectra,  one  gets  a  scaling  factor  which  can  be  directly  convened  to  a  column  density. 
This  method  allows  one  to  eliminate  overlapping  or  contaminating  emissions  from  the 
spectrum  (Cleary,  1986).  Though  not  pursued  in  this  thesis,  the  resulting  fit  can  be 
integrated  over  wavelength  to  get  a  column  emission  rate  for  any  specific  transition.  In 
addition,  the  fit  can  be  converted  to  a  volume  density  through  standard  inversion 
techniques. 

Before  the  synthetic  spectrum  was  fit  to  the  data,  it  was  conditioned  to  model  the 
response  of  the  instrument.  This  involved  convolving  the  synthetic  spectra  with  the 
instrument  slit  function.  As  noted  above,  the  instrument  function  was  found  by  examining 
spectral  lines  of  atomic  mercury.  Figure  5-2  shows  the  effect  of  this  convolution. 

The  11  NO  emissions  identified  and  fit  are  shown  in  Figure  5-3.  The  synthetic 
spectrum  in  this  figure  is  a  summation  of  the  gamma  and  epsilon  bands.  Even  with  the 
low  instrument  resolution,  the  individual  peaks  are  still  distinct  and  relatively 
tmconteiminated  by  other  bands.  A  model  of  the  Nil  2143  doublet  (not  shown)  is  also 
included  in  the  fit  of  the  peak  at  2150  A. 

Not  all  1 1  features  are  fit  at  every  altitude.  At  low  altitudes,  the  stronger  emissions 
saturate  the  detector.  At  higher  altitudes,  the  molecular  nitrogen  Vegard  Kaplan  (VK)  and 
Lyman-Birge-Hopfield  (LBH)  bands  tend  to  obscure  the  weaker  nitric  oxide  emissions. 
This  was  determined  (to  some  extent)  by  fitting  the  data.  These  effects  will  be  discussed 
later. 
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Figure  5-2  Representative  Convolved  Synthetic  Spectrum 
C.  TEMPERATURE  FITTING 

Temperature  was  determined  by  matching  the  shape  of  the  synthetic  7(1,0)  band  to 
that  of  the  data.  There  are  several  factors  which  complicate  this  process.  For  example,  the 
intensity  recorded  by  the  instrument  is  a  line-of-sight  column  emission  rate.  In  the 
spherical  geometry  of  the  atmosphere,  one  would  expect  contributions  to  the  spectrum 
from  higher  altitude,  higher  temperature  regions  along  the  slant  path.  The  effect  is  small 
for  an  exponential  atmosphere,  but  results  in  some  smoothing  of  the  spectmm  from  higher 
temperature  contributions  (see  e.g.,  Bosserman,  1989).  Additionally,  the  wavelength 
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resolution  of  the  instrument  is  such  that  the  rotational  stmcture  of  a  band  is  lost. 
Therefore,  the  sensitivity  of  the  shc^  of  a  given  band  to  changes  in  temperature  is  low. 
Because  of  this,  the  uncertainty  in  determining  a  temperature  fit  was  taken  as  100  K.  The 
net  effect  of  this  uncertainty  is  to  introduce  a  four  percent  systematic  error  into  the 
inferred  column  densities. 

Because  of  the  instrument  resolution,  the  NO  (1,0)  gamma  band  is  blended  with  the 
N*(*S)  doublet  at  2143  A.  For  mid-latitude  observations  at  high  altimdes,  the  NO  doublet 
intensity  is  comparable  to  that  of  the  NO  gamma  (1,0)  band  (Qeary  and  Barth,  1987). 
Therefore,  it  must  be  included  in  the  fit.  The  contribution  from  the  NO  emission  was 
modeled  by  two  delta  functions  at  2143.6  A  and  3139.7  A  with  an  intensity  ratio  of 
100:58  (Buscela  and  Sharp,  1989).  As  with  the  NO  synthetic  spectrum,  the  Nil  model  was 
convolved  with  the  instrument  function. 

The  intensities  of  NO  7(1,0)  and  NH  emissions  were  treated  as  independent 
parameters  in  the  model.  For  a  given  temperature,  these  parameters  were  varied  in  a  grid 
search.  The  synthetic  spectrum  was  fit  to  the  data  so  as  to  minimize  the  chi-squared 
deviations.  The  temperature  used  to  generate  the  NO  synthetic  spectrum  was  then  adjusted 
in  100  Kelvin  steps  until  the  minimum  in  this  third  parameter  was  found. 

A  temperature  profile  was  determined  for  all  altitudes  on  the  up  leg  and  altitudes 
through  145  km  on  the  down  leg.  Below  this,  the  NO  gamma  (1,0)  feature  saturated  the 
detector.  For  the  lower  altitudes,  the  down  leg  temperature  profile  was  extrapolated  by 
comparison  with  the  predictions  of  the  MSlS-83  (Hedin,  1983)  neutral  atmosphere  model. 
Results  of  the  fit  are  shown  in  Figure  5-4.  Temperatures  from  the  data  are  higher  than 
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those  given  by  the  MSIS-83  model,  indicating  higher  solar  activity  than  the  F10.7 
indicates.  As  discussed  below,  a  high  level  of  solar  activity  is  supported  by  the  large 
nitric  oxide  column  densities  observed.  Both  the  up-leg  and  down-leg  temperatures  show 
the  positive  gradient  characteristic  of  the  thermosphere.  The  temperatures  determined  from 
the  fit  were  used  throughout  the  rest  of  the  data  analysis. 


D.  FITTING  THE  DATA 

Having  determined  the  temperature  profile  using  the  y(l,0)  band,  the  remaining  ten 
emission  features  were  fit  with  the  synthetic  spectra.  Each  synthetic  spectrum  was 
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multiplied  by  a  scale  factor  to  match  the  intensity  of  a  specific  emission.  Tl  Js  scale  factor 
is  directly  convertible  to  a  line-of-sight  column  density.  An  average  column  density  for 
a  given  altitude  can  be  calculated  from  the  estimated  column  densities  that  are  determined 
for  the  individual  peaks. 

Not  all  the  emissions  were  suitable  for  the  fitting  process  at  all  aitiiudcs.  as  with 
the  NO  gamma  (1,0)  band,  the  gamma  (0,1),  (0,2),  and  (1,1)  bands  saturated  the  detector 
at  various  altimdes  below  145  km.  In  the  higher  altitude  spectia,  several  of  the  NO  bands 
are  contaminated  by  emissions  from  the  molecular  nitrogen  VK  and  LBH  bands.  The 
wavelength  overlap  of  the  NO  and  Nj  emissions  is  shown  in  Figure  5-5.  Energetic 
electron  impact  is  the  principal  excitation  mechanism  for  the  VK  and  LBH  bands.  The 
photoelectron  flux  is  relatively  low  at  the  lower  altitudes  and  increases  with  increasing 
height  (over  the  altitude  range  considered  here).  Because  of  this,  the  intensity  of  the 
molecular  nitrogen  bands  relative  to  the  nitric  oxide  bands  increases  with  altitude.  The 
overlap  of  the  the  NO  and  Nj  bands  can  distort  the  shape  of  the  NO  emission  spectrum, 
eliminating  distinct  features.  Above  145  km,  only  the  gamma  (1,0),  (2,0),  and  (1,1)  bands 
were  used  in  t;;e  fitting  process. 

The  column  densities  are  plotted  by  v"  progression  for  all  1 1  emissions  in  Figures 
5-6  through  5-9.  The  data  shown  represent  the  down  leg.  Data  firom  185  km  and  195  km 
on  the  up  leg  are  included  to  indicate  trends.  It  is  instructive  to  note  some  of  the  features 
evident.  For  the  gamma  v"  band  progressions  in  Figures  5-6,  5-7,  and  5-8,  there  is  a 
definite  trend  of  Licreasing  inferred  column  density  with  higher  vibrational  states.  The 
relative  intensities  of  emissions  within  a  progression  are  controlled  by  the  Franck-Condon 
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Figure  5-5  Wavelength  Overlap  of  NO  and  Nj  ilmission  Bands 


factors.  Therefore,  the  scale  factors  (or  column  densities)  yield  a  relative  indication  of  the 
correcmess  of  the  Franck-Condon  factors  within  a  progression.  If  the  Franck-Condon 
factors  were  in  the  correct  proportions,  one  would  expect  the  piots  of  the  derived  column 
densities  to  lie  on  top  of  one  another.  Instead,  the  data  shows  a  systematic  trend. 
Recently,  Bosserman,  (1989)  showed  a  similar  trend  in  the  NO  e  bands.  He  proposed 
corrections  for  the  E  band  Franck-Condon  factors  to  bring  his  inferred  column  densities 
into  agreement.  The  Franck-Condon  factors  for  the  e(0,2)  and  e(0,3)  bands  from 
Bosserman  however,  required  little  correction  and  the  present  work  supports  this  (see 
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Figure  5-6  Column  Densities  for  the  y  (0,v")  Progression 
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Figure  5-9  Column  Densities  for  the  E(0,v")  Progression 
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Figure  5-9).  As  a  result,  these  two  bands  were  used  as  a  first  guess  for  the  average 
column  density  profile. 

The  ratios  of  the  y  band  column  densities  to  this  average  column  density  yields  a 
"Franck-Condon  Scale  Factor,"  as  shown  in  Figure  5-10.  These  correction  factors  are 
averaged  over  all  altitudes  and  are  plotted  against  the  band  head  wavelength  for  the 
emissions.  This  plot  shows  that  while  there  is  a  trend  in  the  Frank-Condon  scaling  factors 
within  each  progression,  it  is  not  purely  wavelength  dependent.  Thus,  it  is  unlikely  to  be 
due  to  some  artifact  of  the  sensitivity  calibration  of  the  instrument,  but  rather  is  a  real 
physical  effect. 

A  best  fit  column  density  profile  can  be  obtained  after  correcting  for  these  Franck- 
Condon  scale  factors.  Figure  5-11  shows  the  results  for  the  up-  and  down-leg.  The  error 
bars  in  Figure  5-11  are  plotted  as  three  times  the  uncertainty  for  purposes  of  clarity. 

Figure  5-11  shows  a  difference  of  almost  30%  between  the  up-  and  down-leg 
column  densities  at  155  km.  The  time  between  sets  of  data  was  approximately  five 
minutes,  therefore  it  is  unlikely  that  temporal  variations  in  atmospheric  densities  are  the 
source  of  this  discrepancy.  The  rocket’s  ACS  data  are  stiU  being  analyzed,  but  differences 
in  payload  attitude  may  explain  this  feature  of  the  spectral  data.  Because  of  this,  the  up- 
leg  and  down-leg  spectra  were  treated  separately  in  the  analysis. 

The  column  densities  in  Figure  5-11  are  quite  large  for  mid-latitude  observations 
(see  e.g.,  Cleary  and  Barth,  1986,  McCoy,  1981).  However,  the  Sun  is  approaching  one 
of  the  strongest  solar  maxima  recorded.  The  solar  activity  which  drives  the  production  of 
odd  nitrogen  has  been,  on  the  average,  higher  than  normal.  This  fact  makes  these  large 
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Figure  5-10  Franck-Condon  Scale  Factors 

densities  plausible.  Each  plot  shows  the  characteristic  exponential  altitude  dependence  of 
the  density.  However,  the  upleg  density  falls  off  at  a  slower  rate.  In  the  plot  of  the 
downleg  spectra,  a  peak  near  105  km  is  quite  evident. 

Finally,  Figures  5-12  through  5-24  show  the  best  synthetic  fit  to  the  13  averaged 
spectra.  The  synthetic  fits  have  been  corrected  with  the  Franck-Condon  scale  factor 
adjustments  while  the  wavelength  adjustments  shown  in  Figure  5-1  have  been  added  to 
the  data.  The  data  are  presented  in  the  order  of  their  acquisition;  upleg  spectra  by 
increasing  altitude  and  downleg  spectra  by  decreasing  altitude. 
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Figure  5-11  Average  Column  Etensities 
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Adjusting  the  Franck-Condon  factors  for  the  various  NO  bands  has  clearly  resulted 
in  a  good  fit  to  the  data.  This  is  particularly  evident  at  low  altitudes  where  there  is  less 
contamination  from  other  sources.  At  high  altitudes,  the  contamination  from  the  VK  bands 
at  long  wavelengths  is  quite  evident.  In  particular,  the  shoulders  on  either  side  of  the 
7(0,1)  peak  near  2360  A  are  due  to  the  VK  bands.  The  VK  bands  also  add  to  the  intensity 
of  the  7(0,0),  (0,2),  and  (2,3)  peaks  near  2265  A,  2310  A,  and  2475  A  respectively.  At 
shorter  wavelengths,  between  2000  A  and  2100  A,  one  can  see  the  contributions  from  the 
LBH  bands  affect  the  shape  of  the  7(2,0)  peak  and  the  two  e  bands.  What  is  not 
accounted  for  is  the  discrepancy  between  the  observed  7(1,0)  peak  and  the  synthetic 
spectrum  from  145  km  to  185  km  on  the  up  and  down  legs.  The  synthetic  fit  uses  both 
the  NO  7  and  e  bands  and  the  Nil  emission  at  2143  A.  There  is  some  overlap  of  this 
peak  with  an  LBH  band,  so  this  may  account  for  some  of  the  discrepancy. 

In  the  low  altitude  fits,  contamination  of  the  spectrum  by  Nj  emissions  is  negligible. 
There  are  however,  some  features  which  stand  out.  One  of  these  is  the  detector  saturation 
in  the  7(1,0)  band,  seen  in  Figures  5-22  through  5-24.  This  causes  the  distortion  of  that 
peak.  Additionally,  at  these  low  altitudes,  one  of  the  weaker  NO  emissions,  the  £(1,4) 
band,  can  be  seen  at  2073  A.  This  emission  feature  was  not  fit  with  synthetic  spectra. 
Below  125  km,  the  synthetic  spectrum  appears  to  over  estimate  the  observed  7(2,0) 
emission  at  2050  A.  This  may  be  due  to  self-absorption  in  the  observed  data,  an  effect 
which  has  been  noted  by  other  researchers  (Qeary,  1986).  Last,  there  is  a  divergence 
between  the  observed  intensity  of  the  7(1,1)  band  at  2240  A  and  the  synthetic  fit.  Above 
135  km,  the  7(1,1)  band  synthetic  spectra  provides  the  best  fit  of  all  the  emissions.  Below 
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135  km,  this  changes.  The  synthetic  band  fits  the  data  quite  well  on  the  long  wavelength 
side  of  the  peak,  but  the  peak  intensity  and  short  wavelength  fit  show  a  marked 
difference.  This  suggests  that  the  difference  may  be  caused  by  another  emission,  though 
for  now,  it  remains  unexplained. 
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Figure  5-12  Comparison  of  Data  and  Synthetic  Fit  at  155  km  (Up  Leg) 
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Figure  5-13  Comparison  of  Data  and  Synthetic  Fit  at  165  km  (Up  Leg) 
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Figure  5-14  Comparison  of  Data  and  Synthetic  Fit  at  175  km  (Up  Leg) 
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Figure  5-15  Comparison  of  Data  and  Synthetic  Fit  at  185  km  (Up  Leg) 
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Figure  5-18 


Comparison  of  Data  and  Synthetic  Fit  at  165  km  (Down  Leg) 
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Figure  5-20  Comparison  of  Data  and  Synthetic  Fit  at  145  km  (Down  Leg) 
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Figure  5-22  Comparison  of  Data  and  Synthetic  Fit  at  125  km  (Down  Leg) 
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Figure  5-23  Comparison  of  Data  and  Synthetic  Fit  at  1 15  km  (Down  Leg) 
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VI.  CONCLUSION 


A.  SUMMARY  OF  FINDINGS 

Ultraviolet  emissions  of  the  Earth’s  ionosphere  in  the  wavelength  range  of  2000  A 
to  2500  A  were  analyzed.  Eleven  of  the  most  prominent  NO  gamma  and  epsilon  bands 
in  the  data  were  fit  with  synthetic  spectra.  An  average  column  density  was  found  for 
altitudes  between  100  and  200  km  which  was  larger  than  previous  measurements.  The 
atmospheric  temperature  profile  was  determined  fi’om  the  fit  of  the  (1,0)  gamma  band.  A 
comparison  of  the  temperature  profile  to  that  predicted  by  the  MSIS-83  model  indicated 
a  higher  level  of  solar  activity  than  the  F10.7  would  indicate.  This  conclusion  was 
supported  by  the  exceptionally  high  nitric  oxide  densities  foimd. 

The  derived  column  densities  were  compared  by  v"  progression.  The  data  suggest 
that  the  Franck-Condon  factors  used  in  calculating  the  synthetic  spectra  underestimate 
emissions  of  the  higher  vibrational  states  of  the  gamma  bands.  Scale  adjustments  for  the 
Franck-Condon  factors  were  calculated  and  are  shown  in  Figure  5- 10.  These  scale  factors 
were  applied  to  the  synthetic  spectra  when  determining  the  average  column  density. 

B.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  nitric  oxide  column  densities  found  in  this  work  produce  an  excellent  fit  to  the 
data  for  the  NO  emissions.  There  are  however,  several  areas  of  analysis  which  require 
further  investigation.  For  example,  the  NO  contributions  can  now  be  used  in  the  analysis 


of  other  emission  features  of  the  spectra.  With  the  NO  column  densities  as  a  starting 
point,  Nj  synthetic  spectra  can  be  fit  to  the  data  to  determine  column  densities  for 
molecular  nitrogen.  This  same  technique  can  also  be  applied  to  the  OH  emission  at  2470 
A.  These  steps  are  necessary  to  meet  the  overall  goals  of  the  experiment.  The  low  altitude 
divergence  between  the  observed  y  (1.1)  band  and  its  synthetic  fit  should  also  be 
examined.  As  mentioned  above,  this  difference  may  be  indicative  of  an  overlapping 
emission  from  some  other  source.  The  possibility  of  self  absorption  in  the  low  altitude 
spectra  of  the  y  (2,0)  was  mentioned.  This  effect  could  be  investigated  to  determine  if 
this  is  a  significant  factor  in  this  data.  Unfortimately,  similar  low  altitude  data  for  the  y 
(1,0)  band  is  not  available  from  this  experiment.  Additionally,  completion  of  the  post 
flight  analysis  should  give  some  insight  as  to  the  ^jparent  wavelength  shift  shown  in 
Figure  5-1.  Finally,  future  analyses  of  nitric  oxide  emissions  should  examine  the  observed 
and  calculated  Franck-Condon  factors  for  the  v"  progressions  to  establish  the  trend  noted 
here. 
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